Electron spectra close to a metal-to-insulator transition 
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A high-resolution investigation of the electron spectra close to the metal-to-insulator transition 
in dynamic mean-field theory is presented. An all-numerical, consistent confirmation of a smooth 
transition at zero temperature is provided. In particular, the separation of energy scales is verified. 
Unexpectedly, sharp peaks at the inner Hubbard band edges occur in the metallic regime. They are 
signatures of the important interaction between single-particle excitations and collective modes. 

PACS numbers: 71.30.+h,71.27.+a,71.28.+d,75.40.Gb 



The interplay of electronic degrees of freedom with col- 
lective modes is one of the central issues in current con- 
densed matter physics. It is particularly intriguing when 
the collective modes are formed by the electrons them- 
selves due to strong interactions. Famous examples for 
the complexity of such systems are the high-tcmperaturc 
superconductors (see, e.g., Refs. 1 and 2) and materials 
displaying colossal magnetoresistance (see, e.g., Refs. 3 
and 4). 

We focus on charge and spin degrees of freedom by 
considering a narrow single-band model with nearest- 
neighbor hopping t where the interaction U stems from 
the Coulomb repulsion. For simplicity, we study the half- 
filled case with one electron per site. Thus the minimal 
model is the Hubbard model 



H = -t 



c i;a c j\a 



1/2)^ -1/2) (1) 



where i,j denote sites on a lattice with (i, j) being nearest 
neighbors, a e {j, J,} the spin, c^) a the electron annihila- 
tion (creation), and rii ;t7 their density. 

Leaving aside all effects of long-range order like charge 
or spin density waves the system is metallic for weak 
interaction and insulating for strong interaction. The 
weakly interacting system does not have significant ef- 
fects of collective modes because they are overdamped 
by Landau damping. The strongly interacting system is 
a paramagnetic insulator governed at low energies by the 
collective modes which are the magnetic excitations. The 
charge modes display a large gap of the order of U. Thus 
charge and collective modes are well separated in energy 
so that no significant interplay is to be expected. Hence, 
the regime close to the transition between the metal and 
the insulator is the most likely to be influenced by the 
interplay of single-particle and collective excitations. 

An estimate shows that the two-dimensional supercon- 
ducting cuprates are indeed close to the metal-insulator 
transition. The magnetic coupling J is approximately 
given by 4t 2 /U; empirically, one has J ss i/3 so that 
U/W w 1.5 where W = 8t is the band width in 2D. This 
can be compared to the value of about U /W « 1.2 where 
the paramagnetic insulating phase becomes instable due 
to closing of the charge gap. ! ' 



It is the aim of the present work to provide evi- 
dence that a generic paramagnetic system with values 
of U/W in the range 1 to 1.5 constitutes a highly cor- 
related metal with significant interplay between single- 
particle and collective modes. Since a controlled treat- 
ment of finite-dimensional systems other than the one- 
dimensional chain is not possible we study the Hubbard 
model on the Bethe lattice with infinite coordination 
number z — > oo. 

Scaling the hopping ' t = D/(2y/z) leads to the dy- 
namic mean-field theory (DMFT). 7 > 8 > 9 In DMFT the lat- 
tice problem is mapped onto an effective single impurity 
Anderson model (SIAM) with the self-consistency condi- 
tions that the interaction U, the full local Green's func- 
tion G(io) and the local self-energy S(w) are the same in 
the lattice and in the SIAM. 10 The DMFT as an approx- 
imation to finite-dimensional systems is by now a widely 
employed technique. 11 It allows one to include impor- 
tant correlation effects in the description of real materials 
based on ab initio density-functional theories. Hence the 
quantitative understanding of all features of the DMFT 
solution of the Hamilton operator (1) is mandatory. 

At low, but finite, temperatures T the phase transition 
between metal and insulator is of first order and takes 
place at U C (T) (Refs. 12 and 13) between U C \{T) (insta- 
bility of the insulator against infinitesimal changes in am- 
plitude) and U c2 {T) (instability of the metal). At T = 0, 
the transition has peculiar properties. 14,15 ' 16,1 ' It bears 
features from first-order transitions: a jump in the en- 
tropy and a finite hysteresis between U c \ := U C \{T = 0) 
and U C 2 := U C 2(T = 0). But there are also second- 
order features because U c := U C (T = 0) = U C 2 where 
the quasiparticle weight vanishes continuously and the 
ground-state energy E{U) is differentiablc. The behav- 
ior of E{U) is derived from a projective DMFT, 1S which 
is based on the hypothesis that the energetically high- 
lying spectral features do not change at the transition 
from the metallic to the insulating solution (separation 
of energy scales). 

There are many determinations of U c \ and U C 2- They 
(almost) agree on (2.39 ± 0.02)L> for 1 ' 1,20,21 U cl with the 
outlier of (2.225 ± 0.025)D. 22 For U c2 the values range 
from 2.92D to 3.0D. 16,18,19,21,:23,24 The spectral densities 
p(ui) := — 7r _1 ImG(u;) [density of states (DOS)] display 
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FIG. 1: (Color online) Spectral densities p(uS) deep in the 
metallic regime (upper row) and deep in the insulating regime 
(lower row) in DMFT for the Bethe lattice at T = 0; dashed 
lines: numerical renormalization group (NRG) data (Ref. 32). 



a quasiparticle peak at lj — in the metallic solution 
which is pinned to its noninteracting value po(0). But 
its width decreases on U — ► [/ C 2. 14 ' 16 ' 18 ' 21 The insulating 
solutions display the lower and the upper Hubbard bands 
which merge for U — > E/ C i when the single-particle gap A 
closes. 14 ' 19 ' 21 ' 22 

Recent progress in the numerical calculation of dy- 
namic quantities for quantum impurity models 25 ' 26 ' 2 ' 28 
by dynamic density-matrix renormalization 29 ' 30,3 
(D-DMRG) make calculations possible with well- 
controlled resolution at all energies. Thereby, spectral 
functions and ground-state energies become accessible 
which have so far eluded a quantitative determina- 
tion. With the correction vector method we compute 
p(lj) broadened (convolved) by Lorentzians of width 
r\ 6 [0.01, 0.1] D. The unbroadened p(ui) is retrieved 
by least-bias dcconvolution. 2s It is used to determine 
the continued fraction of the bath function in the next 
iteration of the DMFT self-consistency cycle. 10 For all U 
about 20 iterations were performed till two subsequent 
p(u>) differed less than as 10~ 3 /D everywhere and the 
ground-state energy and the double occupancy differ less 
than 10 _2 %. For the insulator, it is required in addition 
that the static gap, derived from energy differences of 
the finite bath representation, differs less than 1%. 

The D-DMRG is performed with 128 or 256 basis 
states. We use 120, 160, or 240 fcrmionic sites including 
the impurity in the metallic regime. For the insulating 
solutions we used 121 or 161 fermions. An odd num- 
ber of sites implies a pole at lj = in po(u>). This pole 
is split by the interaction. The splitting results from a 
pole in S(w) at lj = 0. Such a solution is insulating. 
Hence an odd number of sites is slightly biased toward 
an insulator. Vice versa, an even number of sites leads 
to ImE(0) = implying a small bias toward the metallic 
solution. The relative bias is estimated by the inverse 
number of sites: (4-8) xl0~ 3 . In odd chains, we observe 
two almost degenerate ground-states (spin f or J, at the 
interacting site) which must both be considered. Other- 
wise a spurious magnetic moment is generated. 




co/D 

FIG. 2: (Color online) Spectral densities p(u>) of the metallic 
(solid) and the insulating (dashed) solutions between U c i and 
U c2 . 



In Figs. 1 and 2, our results for metallic and insulating 
p{u>) are shown. In the metallic solutions, the narrowing 
of the quasiparticle band around lj = is clearly visible. 
From U ~ D on, the DOS displays side features which 
develop into the lower and upper Hubbard bands. At 
U ~ 2D the Hubbard bands are well separated from the 
quasiparticle peak at lj = by a precursor of the gap 
A in the insulator: a pseudogap. The comparison with 
the NRG data from Ref. 1 ( i shows good agreement in the 
quasiparticle peak but deviations in the Hubbard bands. 
There the DMRG data are much sharper and do not have 
significant tails at higher energies. This difference stems 
from the broadening proportional to the frequency which 
is inherent to the NRG algorithm. 25,26 

The insulating solutions display the lower and the up- 
per Hubbard bands clearly. They agree excellently with 
the perturbative result 33 (not shown) for U ^ 3D. At 
U = U c i = (2.38 ± 0.02)D both bands touch each other. 
No upturn in p(uj) as in Ref. 22 is found when we con- 
sider the deconvolved p(ui) for all lj. An upturn occurs 
only if the static gap is used. But such a procedure did 
not lead to stable self-consistent solutions. 

In Fig. 3 the quasiparticle weight Z in the metal and 
the single-particle gap A in the insulator are shown. The 
weight Z = [1 — 9 w ReS(0)] _1 is found from fitting the 
derivative of the Dyson equation G(lu) = Gq[uj — £((*;)] 
implying Z~ x = D 2 d UJ G(0)/2 where Gq{u) is the bare 
local Green's function of the lattice. The gap A is found 
from a fit proportional to \/oj — A to the sharp uprise 
of the DOS at the inner band edges. The DMRG data 
agree excellently with the perturbative results where the 
respective perturbation holds. The comparison to the 
NRG data shows that the broader high energy features 
reduce Z to some extent so that the NRG weight stays 
below the DMRG data. 

From the power-law fit A = (U— J7 c i) ? [ai+a 2 (U— £7 c i)] 
shown in Fig. 3 we find U c i = (2.38 ± 0.02)£> in perfect 
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FIG. 3: (Color online) Dotted area: region of two solutions. 
Left curves: metallic quasiparticle weight Z; line with circles: 
interpolated DMRG, line with pluses: NRG 11 '; dashed line: 
perturbation up to U 4 (Ref. 25). Right curves: insulating 
gap A or pseudogap in the metal (line with diamonds); line 
with squares: DMRG; dashed line: perturbation up to 1/U 2 
(Ref. 33). Inset: weight S of the peaks at the inner Hubbard 
band edges. 



agreement with most of the previous results. 19,20 ' 21 The 
exponent is found to be £ = 0.72 ± 0.05. The value for 
U C 2 = (3. 07 ±0.1)1) is determined reliably from a second- 
order extrapolation shown in Fig. 3. The value of U C 2 
agrees well with the previous results. 16 ' 18 ' 19,21 ' 23 ' 24 We 
attribute the small deviation to the enhanced accuracy 
of our D-DMRG algorithm. 

Figure 2 nicely shows the evolution of the metallic and 
the insulating p{w) between U c \ and U C 2- It represents 
the first all-numerical high-precision confirmation of the 
hypothesis of the separation of energy scales. Clearly, the 
metallic Hubbard bands at higher energies approach the 
insulating ones for U — > U C 2- Differences remain at the 
inner band edges and around u> = as long as U < U C 2- 
The frequency of the sharp uprise in p{u>) at the inner 
edges in the metallic solutions defines the pseudogap. Its 
evolution (line with diamonds in Fig. 3) shows that it 
tends to the insulating gap at U = U C 2- This corrobo- 
rates strongly that the metallic p(oj) approaches the in- 
sulating p{ui) for U — » U C 2- Remarkably, the metal and 
the insulator have the same single-particle correlations 
at U = U C 2, which is again a feature of a second-order 
phase transition. 

From equations of motion"^ one obtains the ground- 
state energy per site 



E(U) 
L 



E 



.t 



[c, (T ,H])-U[d{U) 



(2) 



where d(U) is the double occupancy (rn^n^i) . All the 
local expectation values are computed directly at the 
interaction site of the SIAM at self-consistency. Then 
they equal the local expectation values on the lattice. 
In Fig. 4, E(U) and d(U) are depicted. The equality of 
d(U) in the SIAM [data points and interpolations (solid 
lines)] and in the lattice, obtained from d(U) = djjE+1/4 




FIG. 4: (Color online) Upper panel: ground-state energy 
E(U) in the insulator (squares) and in the metal (circles). 
Lower panel: corresponding double occupancy d(U). At U c = 
(3.07 ± 0.04)D the double occupancies intersect. The dashed 
lines depict the derivative of the interpolations for E{U). 



(dashed lines) , is a very sensitive check of the data qual- 
ity. It is perfectly fulfilled in the insulator (difference 
$2x 10~ 4 ). The comparison of energies after integrat- 
ing d(U) for U € [2.4, 4.0}D as in Ref. 22 shows agreement 
within 8 x 10 _6 D. In the metal, our check for d(U) works 
very well, but not perfectly (difference ~ 3 x 10~ 3 ). The 
remaining deviation is a finite-size effect of the bath be- 
cause the deviation decreases on increasing the number 
of fcrmionic sites. 

Assuming a differcntiable behavior of E(U) we deduce 
the critical interaction U c = (3.07 ± 0.04)£> from the in- 
tersection of the double occupancies. The corresponding 
metallic and insulating energies differ less than 10" 4 D. 
This justifies the assumption of differentiability a pos- 
teriori. The agreement of U c and U C 2 in our numerical 
treatment proves the consistency of our data. Thereby, 
the previously proposed scenario 14 ' 15 ' 18 for the metal- 
insulator transition is numerically confirmed. 

An unexpected feature in our metallic solution are the 
sharp peaks at the inner edges of the Hubbard bands 
(see Fig. 2 and Fig. 1 for U = 2D). The only previ- 
ous evidence for similar features were weak shoulders in 
NRG 16 and quantum Monte Carlo 1 ' 5 data. Based on the 
following arguments we interpret it as the signature of a 
collective mode. 

Since the sharp peaks occur only in the metal the heavy 
quasiparticles must be involved. But the peaks are lo- 
cated at relatively high energies compared to the quasi- 
particle band. So a binding or antibinding phenomenon 
with something else must take place. Since the peaks 
are very sharp (to numerical accuracy for U 2.61?) we 
conclude that an immobile, local mode is involved. Spin 
excitations like spin waves are the best-known collective 
excitations in Hubbard systems. In infinite dimensions, 
they are indeed dispersionless. 34 In absence of any mag- 
netic order they are located at zero energy in the insu- 
lator because the Hartree term vanishes. ' By continuity, 
we deduce that their energy is very low also in the metal 
with an estimated upper bound of 0.2D deduced from the 
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energy of a spin wave in the long-range ordered Neel state 
which is zJ/2 ps D 2 /(2U) rts 0.2D at U = 2.5D. Hence, 
we conclude from the energy of the sharp features that 
it signals the antibound state or resonance of a heavy 
quasiparticle with a collective spin excitation. We refer 
to this state as the antipolaron. 

In the inset of Fig. 3 we plot the weight S of the an- 
tipolaron peak as a function of Z. The error bars re- 
sult from the numerical difficulty to resolve this sharp 
feature and from the analytical difficulty to separate it 
from the background of the Hubbard bands. For the 
separation, we fitted an approximate square root onset 
starting at the pseudogap, multiplied with a quadratic 
polynomial, to the Hubbard bands. The excess weight 
S is attributed to the antipolaron peak. From the inset 
in Fig. 3 we conclude that S vanishes linearly with Z, 
rather than quadratically or cubically. Thus the antipo- 
laron peak vanishes proportional to the matrix clement 
of a single quasiparticle. If S had vanished quadratically 
(cubically) one would have concluded that two (three) 
quasiparticles were involved. The linear dependence of 
S on Z corroborates the scenario of an antibound state 
formed from one quasiparticle and one collective mode. 
Surely, further investigations are called for. 

We point out that the antibinding between the heavy 
quasiparticles and the collective magnetic modes suggests 
an interesting answer to the so far open question why the 
metal forms Hubbard bands and eventually becomes in- 



sulating. The weight close to the Fermi level is pushed 
away to higher energies by a strong repulsive interac- 
tion between the low-lying quasiparticles and magnetic 
modes. Note that this scenario is also applicable in finite 
dimensions if the finite dispersion of the collective modes 
is accounted for. 

In summary, we provided the all-numerical con- 
firmation of the metal-insulator transition proposed 
earlier 14 ' 15 ' 18 on the basis of the hypothesis of the sepa- 
ration of energy scales. We found that the spectral den- 
sity of the metal approaches the one of the insulator for 
U -» U c2 . The critical interactions U cl = (2.38 ± 0.02)D 
and U C 2 = (3.07 ± 0.1)D were found; the latter value co- 
incides with U c = (3.07 ± 0.04)£> where the ground-state 
energies E(U) intersect differentiably. Hence, a consis- 
tent picture of a diffcrentiable E(U) emerged. 

The unprecedented resolution of the spectral densities 
enabled us to detect sharp peaks at the inner Hubbard 
band edges. We interpreted them as antibound states 
or resonances formed from a heavy quasiparticle and a 
collective magnetic mode. The occurrence of such signa- 
tures in electron spectra opens up an interesting route to 
investigate the interplay of single-particle and collective 
excitations in photoelectron spectroscopy. 
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sions and the DFG for financial support in SFB 608. 
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